Gauss sum factorization with cold atoms 
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We report the first implementation of a Gauss sum factorization algoritlim by an internal state 
Ramsey interferometer using cold atoms. A sequence of appropriately designed light pulses interacts 
with an ensemble of cold rubidium atoms. The final population in the involved atomic levels 
determines a Gauss sum. With this technique we factor the number 7V=263193. 



The Shor algorithm [l[ to factor numbers and its NMR- 
implemcntation [21 have propelled the field of quan- 
tum cornputation [3|. Recently a different factorization 
scheme [J] taking advantage of the periodicity proper- 
ties of Gauss sums Q has been proposed and veri- 
fied by two NMR-experiments 7|, |8| and one experiment 
based on short laser pulses [9|]. In the present paper 
we report the first implementation of Gauss sum factor- 
ization based on matter-wave interfcrometry [lo| with 
cold rubidium atoms and use it to find the factors of 
N = 263193. 

Our method rests on the observation that the trun- 
cated Gauss sum 
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consisting of M + 1 terms yields unity if the integer / is 
a factor of A^. In contrast, for integer non- factors I de- 
structive interference leads to a small value of C^j^\l). 
Thus the algorithm checks if a given integer / is a factor 
of N or not. This testing is performed by a quantum sys- 
tem and its time evolution determines if a trial factor is a 
factor or not. Since in the worst case wc have to test all 
prime numbers up to \fN the prime number theorem Q 
predicts the upper bound ^/N /\ogN for the number of 
trials. As a consequence our method in the present form 
scales exponentially which is not surprising since it does 
not involve entanglement yet and relies solely on interfer- 
ence. In this respect, our technique is very much in line 
with the recent critical discussion llj of the connection 
between quantum mechanics and factorizing. 

Our procedure for implementing Gauss sums is rem- 
iniscent of the method used in Ref. 0. Common to 
both techniques is a sequence of pulses which imprints 
on a two-level quantum system a sequence of well-defined 
phases. For appropriately chosen pulses the excitation 
probability takes the form of a Gauss sum. 

Despite these similarities our technique is different in 
many aspects: (i) In contrast to the NMR-approaches 0, 
Q we do not use liquids but an ensemble of cold atoms. 
The spectacular control of internal as well as external 
atomic degrees of freedom with the help of lasers, opens 
up a new avenue towards factorization. Indeed, rubidium 



FIG. 1: Implementation of the Gauss sum factorization al- 
gorithm using cold rubidium atoms. They are launched by 
a double MOT system and prepared in the region P by an 
appropriate pulse sequence in the atomic ground state. In 
the factorization zone F the atoms interact with a sequence 
of pulses driving a hyperfine transition. We start with a 7r/2- 
pulse followed by a sequence of m + 1 7r-pulses with the phase 
4>k{l) and conclude by another 7r/2-pulse. The phase during 
the 7r/2-pulses is = (jyj = —90°. A fluorescence detec- 
tion in the region D measures the populations in both states 
which determines the interference signal c,„(i). We record Cm 
for < m < M and the sum over Cm yields the Gauss sum 
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atoms feature two long-living hyperfine ground states 
which can be coherently manipulated by a two-photon 
Raman-transition [l2|. Moreover, the use of cold atoms 
provides us with long interaction times and the possibility 
of a large number of pulses, {ii) In our approach there are 
no projection measurements between the pulses Q- The 
measurement takes place after completing the sequence 
of pulses. {Hi) The Gauss sum results from the addition 
of the measured signals obtained from an increasing num- 
ber of pulses. As a consequence our resources in pulses 
scale quadratically whereas in the NMR-approach they 
only scale linearly. Despite this unfavorable scaling of 
the present technique we maintain that cold atoms offer 
several advantages. Most importantly, they have already 
proven to be ideal objects to be entangled. In particu- 
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lar, two-qubit quantum gates relying on cold atoms in 
standing light waves have been realized. Therefore, 
we consider our experiment as a stepping stone towards 
more complex arrangements including entangled quan- 
tum systems. 

Our experimental setup shown in Fig. [T] is part of an 
atom interferometer [14| designed to measure rotations 
with a high precision. The source of the cold Rb-atoms 
IS a double stage MOT [H, A two-dimensional trap 
creates an atomic beam with a flux of 5 x 10^ at/s for 
the efficient loading of a subsequent 3D-M0T. Approxi- 
mately Nat = 10^ trapped atoms are launched in a mov- 
ing molasses, similar to atomic fountains [l7l | and arc 
transferred into the atomic ground state 15^5*1/2, F — 
1, mi? = > using a multi-stage preparation. 

In order to implement the Gauss sum, we use the pulse 
sequence illustrated in Fig. [TJ An initial 7r/2-pulse with 
phase (pi = —90° prepares a coherent superposition of 
ground and excited state. The latter corresponds to the 
hyperfine ground state 15^5*1/2, F = 2, = >. This 
superposition is most sensitive to the phases of 7r-pulses. 
Indeed, after a time T we apply the factorization se- 
quence consisting of m -|- 1 7r-pulses each separated by 
the time T. The fc-th pulse has the phase 0] 







for fc = 



ak{N)/l for 1 < k 



(2) 



with ak{N) = {—l)^TTN{2k — 1) which induces a coher- 
ent evolution of the atomic ensemble. We conclude at 
the time (m -I- 1)T the factorization sequence by a 7r/2- 
pulse with phase 0/ = —90° to convert the phase evolu- 
tion into a population difference between the two atomic 
states which is measured by a state-selective fluorescence 
detection [lit . 

The populations in the two states are governed 0, 0| 
by the interference signal 



c„i(/) = cos 2TTm 
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which assumes values between -1 and -1-1 corresponding 
to all atoms in the ground or excited state, respectively. 
When we repeat the pulse sequence for increasing m with 
< m < M and add the interference signals we arrive 
after normalization at the total signal (l) determined 
by Eq. ©. 

The multi-pulse excitation is executed by two digitally 
phase-locked Raman-lasers [3| which drive the hyperfine 
transition at approximately 6.834 GHz. The two beams 
are co-propagating through the factorization zone per- 
pendicular to the trajectories of the atoms. In contrast 
to the velocity-sensitive excitation in inertial atomic sen- 
sors j2ol | , in the present velocity- insensitive configuration 
we can neglect phase contributions from inertial forces. 

A low phase noise oscillator serves as the reference for 
the phase locking setup. The phase 4>k{l) determined in 



advance by a computer for each trial factor I is adjusted 
electronically by a synthesizer. The stringent require- 
ments on the phase control in atom interferometry make 
the realization of the factorization experiment possible. 
The total phase error of the Raman-laser system is ap- 
proximately 1 mrad. The length of a 7r-pulse is approx- 
imately 23 fj,s while the time T between two 7r-pulses is 
100 fis which is sufficient for the electronic control loop 
to adjust the required phase in the laser system. 

We now turn to the discussion of the results of our 
factorization experiment exemplified by the number N = 
263193 = 3x 7x83x 151. Figure[2brings out most clearly 
that the dependence of the interference signal c„i(l) on m 
is dramatically different for a factor of N such as I = 151 
and a non-factor such as I = 150. Since for a factor Cmil) 
is approximately constant as a function of m the con- 
structive addition of the terms in the definition Eq. ([1]) 
of the Gauss sum leads to a total signal C^^^\l) close to 
unity. On the other hand a non-factor produces an os- 
cillatory function Cm(l) which takes on values between -1 
and -|-1. Thus the total signal C^^\l) for a non-factor is 
rather small. 

Figure [2] also shows a decay of the interference sig- 
nal with increasing m. This behavior originates from 
the fact, that the atomic cloud on its trajectory sees 
a Gaussian intensity distribution of the Raman-laser 
beams which translates into a Gaussian distribution of 
the lengths of the tt/2- and 7r-pulses. This leads to a re- 
duction of the transition probability and to a decrease of 




FIG. 2: Measured interference signals Cmil) as a function of 
the summation index m of the Gauss sum Eq. (HI determined 
by the number m -f 1 of pulses in the factorization sequence 
for N = 263193 = 3 x 7 x 83 x 151. Whereas Cm{l) is approx- 
imately constant for the factor 1=151 leading to the value 
C'J^^ (151)=0.69 the corresponding signal for the non-factor 

(14), 
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1=150 oscillates creating the small value CjY*^(150)= -0.02 



Here we have used the maximum number 71/ -|- 1 
pulses. 
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FIG. 3: Average of the interference signals Cm{l) correspond- 
ing to three factors of = 263193 as a function of m with 
and without a parabolic adaption of the pulse length. 

the interference signal during the multi-pulse sequence. 

We have compensated this disturbance by appropri- 
ately adjusting the pulse length over the total interaction 
region. For this purpose we have measured the lengths 
of the TT-pulses for atoms in the center and in the edge 
of the interaction region which is 20 fj,s and 26 fj,s, re- 
spectively. Since we can approximate a Gaussian in the 
neighborhood of its center by a polynomial of second de- 
gree we can connect the two measured pulse lengths by 
a parabola which provides us with an approximation of 
the pulse lengths over the whole interferometer region. 

In Fig. [3] we compare the interference signals with and 
without such a pulse length adaption in their dependence 
on m. For this purpose we concentrate on factors only. 
In order to obtain a smooth curve we take the average of 
three such curves each corresponding to a different fac- 
tor. This procedure clearly shows that with the adaption 
the decay is slowed down and consequently we can apply 
more pulses. 

The remaining decay of c„j(/) results from the inter- 
action of the atomic cloud which has currently a diam- 
eter of about 5 mm in the interaction region with the 
Gaussian intensity distributed Raman-lasers leading to a 
reduction of the transition probability during each single 
pulse. Thus lower temperatures of the atoms would slow 
down the decay of Cm{l)- 

An upper limit for the number M -|- 1 of possible pulses 
in the factorization sequence is experimentally set by the 
total interaction time [M -\- 2)T which results from the 
non- vanishing velocity Vat ~ 4.4 m/s of the atoms prop- 
agating through the interaction region given by the diam- 
eter d ^ 30 mm of the Raman-laser beams. Thus we can 
apply up to 20 pulses without any excessively disturbing 
effects. 

Figure 3] demonstrates the successful implementation 
of a Gauss sum factorization algorithm by an internal 



FIG. 4: Experimental Gauss sum factorization based on cold 
atoms and exemplified by the number A'^ = 263193 = 3 x 
7 X 83 X 151 using the maximum number M + 1=15 of pulses. 
Here we display the absolute value of the sum C^^' (l) over the 
interference signals Cm(Z), that is the measured Gauss sum as 
a function of the trial factors 1 < / < 200. Dominant peaks 
correspond to factors of N. The background is produced by 
the non-factors. 

state Ramsey interferometer using cold atoms. Here we 
display the absolute value of the sum C\^'\l) of the inter- 
ference signals Cm(0 for analyzed integer trial factors 
/ for the number A^ = 263193 = 3 x 7 x 83 x 151 recorded 
with the maximum number M -|- 1 = 15 of pulses. The 
corresponding factors of A^ stand out clearly from the 
background created by the non-factors. We also obtain 
factor-like signals for products of factors as exemplified 
by 1=21, in complete agreement with the Gauss sum, 
Eq. HI. 
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FIG. 5: Measured contrast V of the factorization pattern for 
A*'=263193 as a function of M. Here M -I- 1 is the maximum 
number of factorization pulses. Already six pulses (Af = 5) 
yield a contrast larger than 60%. 
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Figure [5] displays the experimentally obtained contrast 
V of a factorization pattern such as the one in Fig. [4] as 
a function of the number M. Following Ref. 0] we have 
defined V as the ratio between the difference and the sum 
of the measured averaged absolute values of the Gauss 
sum Cj^f^^ (l) at factors and non-factors. We find that with 
just six pulses corresponding to Af = 5 the number N = 
263193 can be factored with a contrast larger than 60%, 
in agreement with the theoretical prediction of Ref. 0] . 

We have demonstrated the first implementation of 
a Gauss sum algorithm to factor numbers based on 
cold atoms and have successfully factorized the number 
A^=263193. Our results are comparable to those of the 
NMR-expcriments 0, S|- However, the use of cold atoms 
not only represents an alternative approach but also al- 
lows us to envision several extensions: (z) Quadratic 
phases can be generated without prior calculation by a 
computer by a linear sweep of an external magnetic or 
electric field during the multi-pulse sequence (?«) The 
preparation of the atomic ensemble in an optical lattice 
opens up the possibility of applying an almost arbitrarily 
large number of pulses. (Hi) A Gauss sum factorization 
scheme involving entanglement could rely on a large num- 
ber of entangled atoms each one located in the minima 
of an optical lattice [21| providing us with a massive par- 
allelism. The experiment reported in the present paper 
is the first step in these directions. 
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